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Three recently suggested random matrix ensembles (RME) 
are linked together to represent a class of RME with multifrac- 
tal eigenfunction statistics. The generic form of the two-level 
correlation function for the case of weak and extremely strong 
multifractality is suggested. 
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Random matrix ensembles turn out to be a natural 
and convenient language to formulate generic statistical 
properties of energy levels and transmission matrix ele- 
ments in complex quantum systems. Gaussian random 
matrix ensembles, first introduced by Wigner and Dyson 
' for describing the spectrum of complex nuclei, be- 
came very popular in solid state physics as one of the 
main theoretical tools to study mesoscopic fluctuations 3 
in small disordered electronic systems. The success of the 
random matrix theory (RMT) 2 in mesoscopic physics is 
due to its extension to the problem of electronic trans- 
port based on the Landauer-Biittiker formula 4 and the 
statistical theory of transmission eigenvalues 5 . Another 
field where the RMT is exploited very intensively is the 
problem of the semiclassical approximation in quantum 
systems whose classical counterpart is chaotic 6 . It turns 
out 6 that the energy level statistics in true chaotic sys- 
tems is described by the RMT, in contrast to that in the 
integrable systems where in most cases it is close to the 
Poisson statistics. 

Apparently the nature of the energy level statistics 
is related to the structure of eigenfunctions, and more 
precisely, to the overlapping between different eigenfunc- 
tions. This is well illustrated by spectral statistics in a 
system of non-interacting electrons in a random potential 
which exhibits the Anderson metal-insulator transition 
with increasing disorder. At small disorder the electron 
wave functions are extended and essentially structureless. 
They overlap very well with each other resulting in energy 
level repulsion characteristic of the Wigner-Dyson statis- 
tics. On the other hand in the localized phase electrons 
are typically localized at different points of the sample, 
and in the thermodynamic limit where the system size 
L — > oo they "do not talk to each other". In this case 
there is no correlation between eigenvalues, and the en- 
ergy levels follow the Poisson statistics. 

The energy level statistics in the critical region near 
the Anderson transition turns out to be universal and dif- 
ferent from both Wigner-Dyson statistics and the Pois- 



son statistics 7,8 . A remarkable feature of the critical 
level statistics is that the level number variance £2 (N) = 
{(SN) 2 ) = \N is asymptotically linear in the mean num- 
ber of levels N 3> 1 in the energy window. Such a quasi- 
Poisson behavior was first predicted in Ref. [9], Later 
the existence of the linear term in T,2(N) was questioned 
8 , since for this term to appear the normalization sum 
rule should be violated. It has been shown recently 10 
that the new qualitative feature responsible for the vi- 
olation of the sum rule and the existence of the finite 
"level compressibility" \ is the multifractality of critical 
wave functions ll > 12 . 

The notion of multifractality is two-fold. The first (and 
most widely accepted) property of multifractality is re- 
lated with the space structure of a single wave function 
*S> n (r) . It is defined through the moments of inverse par- 
ticipation ratio 11 : 

^ = E(i vI '«( r )i 2p ) (xL ^ p(p ~ 1) ' (!) 

r 

where L is the system size, d is the dimensionality of 
space; p > 1 is an integer. The set of exponents D p < 
d characterize the fractal dimensionality of the cluster 
where |^ r „(r)| is larger than a certain value that increases 
with increasing p. 

The second, far less appreciated property of multifrac- 
tality, is related to the overlapping of different wave func- 
tions with energies E n and E m . The main effect of mul- 
tifractality on spectral statistics is given by the simple 
overlapping of local densities (p — 2). The correspond- 
ing fractal dimensionality D2 is the most important crit- 
ical exponent. For \E n — E m \ 3> A (A = l/pL d , where 
p = (p(E)) is the mean density of states) the form of the 
local density correlation function has been suggested and 
confirmed numerically in Ref. [12]: 

(|* n (r)| 2 |v]/ m (r)| 2 ) « \E n - E m \-^\ (2) 

A remarkable feature of multifractality is that the local 
density correlation function decreases very slowly with 
increasing \E n — E m \ so that two fractal wave functions, 
however sparse they are, should still overlap strongly 13 . 

One of the consequences 10 of Eq.(||) is the anoma- 
lous Poisson- like term in the level number variance ^(N) 
which is characterized by the level compressibility \'- 

d-D 2 

X = ^2d— (3) 
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It is immediately seen from Eq.Q that the critical level 
compressibility never reaches the Poisson limit x = 1. 
For an infinitely sparse fractal d — D 2 — > d is maximal, 
yet x is equal to 1/2 and not to 1. This is because even 
the infinitely sparse critical fractal eigcnfunctions overlap 
strongly, in contrast to two localized states 10 < 13 . 

One may assume that the universal critical level statis- 
tics which is described by a set of critical exponents D p , 
applies to a wider class of physical systems and it is in 
fact generic for an intermediate situation between chaos 
and intcgrability. An example of such a system has been 
proposed recently 14 . It turns out that the Coulomb 
impurity inside an integrable square billiard leads to a 
drastic reconstruction of eigenstates however small is the 
strength of the potential. In such a " Coulomb billiard" 
all eigenfunctions in the momentum representation ex- 
hibit multifractality. 

It is therefore natural to look for a RME with mul- 
tifractal eigenvector and eigenvalue statistics similar to 
that at the mobility edge in disordered electronic sys- 
tems. Such a RMT would provide a description of generic 
features of the critical level statistics. 

One such ensemble is suggested in Ref.[15]. It is the 
Gaussian ensemble of M x M hermitian matrices H with 
independent random entries (i > j) defined by: 



<ff«>=0, ((H* i ) 2 )=l3- 1 [l + \i-3\ 3 /B i 



(4) 



where HQ are real random components (fj, = 1 for i = j, 
(i = l, ■■■f3 for i > j); f3 = 1,2,4 for Dyson's orthogonal, 
unitary and symplectic ensembles, and B is a parameter. 
For B»l this RME can be mapped onto a nonlinear 
supersymmetric sigma- model 15 . The case B <C 1 cor- 
responds to 2D Coulomb billiard considered in Ref. [14]. 
The presence of multifractality, Eq.(|l]), and Eq.(||) has 
been proved 14 ' 15 for this RME. In what follows we will 
use this RME (RME-I) as a reference point. 

There are two more RMEs 16,17 which were suggested 
recently as possible candidates to describe the critical 
level statistics. However, their definitions are so dras- 
tically different that they were considered as two alter- 
native options albeit the two-level correlation functions 
(TLCF) R{e,s) = (p(e)p(e + s)) c in the proper regimes 
are identical for both RMEs. It was first noted in Ref. [18] 
that since the energy level statistics is a "fingerprint" 
of the statistics of eigenfunctions, the latter in the cor- 
responding regimes of these two models should also be 
similar. 

The first quantitative link between the predictions of 
RME equivalent to that studied in Ref. [16] (RME-II) and 
numerics on the 3D Anderson model at the mobility edge 
has been done in Ref. [19]. Surprisingly enough it was 
possible to fit very well the numerics for the critical level 
spacing distribution P(s) in 3D Anderson model by a 
proper choice of only one parameter in RME-II. More- 
over, the level compressibility x m the RME-II for this 
particular choice of parameter turned out to be very close 
to that found numerically for the 3D Anderson model. 



In this Letter we argue that RME-II and a certain 
critical regime in RME-III studied in Ref. [17] are equiv- 
alent to RME-I and thus possess the multifractality. Al- 
together they form a new class of RME which describes 
certain remarkable features of critical level statistics. 

We start with the definitions of RMEs studied in 
Ref.[16,17]. The probability density P{H) for a M x M 
random Hamiltonian H from RME-II is given by: 



Pn(H) oc exp[-0TrV(H)], 



(5) 



where the "confinement potential" V(H) grows ex- 
tremely slowly with H: 



V(H) = -\n 2 H, #>1. 
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(6) 



This is crucial for the universality of the eigenvalue statis- 
tics in the limit M — > oo 20 . For V(H) growing slower 
than H the full universality is no longer present, 21:22 
and the eigenvalue statistics may, and does differ from 
the Wigner-Dyson statistics 16 . Another important fea- 
ture of Eq.(|5|) is that the distribution function Pu(H) is 
invariant under the rotation of basis (unitary invariance): 



Pa(H) = Pu{UHtf). 



(7) 



In contrast to Eq.(||), the distribution function for 
RME studied in Ref. [17] is Gaussian. However, the uni- 
tary invariance is broken by a fixed unitary matrix SI: 



Piu{H) <x e-* 



TrH 2 -f &Tr{[n,lf][n,ff] + } 



(8) 



The properties of this RME depend on the choice of fl. 
For the reasons discussed below we consider as RME-III 
the RME defined by Eq.(|) with Q = diag(e l6 i), where 
6j = 2-njjM. The relevant critical regime for this RME 
corresponds to the symmetry breaking field b that scales 
with M — > co as b — h 2 M 2 , where h is a parameter. 

Now it is clear why the RMEs given by Eqs.(||),(||) 
seem so drastically different. The lack of unitary invari- 
ance of Pin(H) means a preferential basis. The existence 
of such a basis implies a certain structure of eigenfunc- 
tions (in this basis) which should lead to spectral statis- 
tics different from the Wigner-Dyson one. However, it 
seems there is no way to get any structure of eigenfunc- 
tions in the unitary-invariant RME-II. It follows imme- 
diately from Eq.(Q) that the distribution function Pn(H) 
depends only on E n , and the statistics of eigenfunctions 
in RME-II is trivial and the same as for standard Gaus- 
sian ensembles 2 . Then the physical picture that the 
spectral statistics is related to the statistics of overlap- 
ping eigenfunctions seems to leave only one single option: 
the Wigner-Dyson energy level statistics in RME-II. 

Nonetheless TLCF R(e,s) = S(s) + Y 2 (e,s) proves to 
be identical in these RMEs and after unfolding 23 it takes 
the form 16 < 17 : 



Y 2 (s,s) 



it r\ 



sin 2 \ns] 



4 sinh 2 [s7r 2 ?7/2]' 



09 = 2), 



(9) 



2 



where r) = j/ir 2 < 1 or r\ = h < 1 for RME-II and RME- 
III, respectively, and e ^> \s\. Eq.(^) coincides with the 
density correlation function for a free electron gas at a 
finite temperature r\£p and differs from the RMT result. 

The way out from this contradiction is suggested in 
Ref. [18] where it was conjectured that the unitary invari- 
ance is broken in RME-II spontaneously. This means that 
the statistics of eigenfunctions in this ensemble should be 
calculated after an infinitesimal symmetry-breaking term 
similar to that in Eq,(||) is added. Then the identical 
TLCF in RME-II and RME-III should be considered as 
an evidence that the proper procedure should result in 
similar eigenfunction statistics in RME-II and RME-III. 

The progress 17 in studying the level statistics in RME- 
III that lead to Eq. (^) is due to averaging over the unitary 
group fi. The level statistics depend on the configuration 
of eigenvalues e l8j of il. The main contribution to the 
average is made by f2 with the most homogeneous con- 
figurations of 9i, the property known as an eigenvalue 
repulsion 2,17 Therefore, one may expect that the spec- 
tral statistics obtained after such an averaging is close 
to that corresponding to a single unitary matrix f2 with 
eigenvalues ftj = exp[(27ri/M)j] (RME-III). As a matter 
of fact for h <C 1 it turns out to be the same. 

In order to show that we note that in the limit M — > oo 
Eq.(|) leads to: 
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TvP" 



1-J\ 



(10) 



If b/M 2 — + 0, then we have a standard Gaussian ensemble 
1,2 and the Wigner-Dyson statistics. In the opposite case 
b/M 2 — > oo, we have an ensemble of random diagonal 
matrices and the Poisson statistics. In the critical case 
considered here for b — h 2 AI 2 , the behavior of ((H^) 2 ) 
is the same as in Eq.(4) defining RME-I. We conclude 
that the M -> oo limits of RME-I and RME-III coincide. 
Then TLCF for RME-III and RME-I must be identical. 

Fortunately, the latter can be calculated directly. 
TLCF can be expressed 24 in terms of the spectral deter- 
minant P{s) as follows: 



Y 2 (s) = 



where 



1 



1 d 2 



In Pis) 

2 C 2 An 2 ds 2 {> 



C0S(27TS) 



2tt 2 s 



p( S ) = n 1 1 



(11) 

(12) 



and e n is a spectrum of the quasi-diffusion modes. The 
latter can be found from the mapping 15 onto the nonlin- 
ear sigma-model (B ^> 1) as follows: e„ = 4B\n\, where 
for the periodic boundary conditions n = ±1, ±2, .... 
Making use of the identity x _1 sinhx = Il^Li(l + 
x 2 /n 2 n 2 ) we immediately arrive at Eq.(^|) with r\ — 
l/(2nB). Using the results of Ref. [15] one can express 



the multifractality exponent Z?2 in terms of B ^> 1. For 
(3 = 2 it appears to be D2 = 1 — 1/(2ttB) which helps to 
identify the parameter rj in Eq.(||) as i] = 1 — Lh. Thus 
all three ensembles share the same TLCF, Eq.([|) which 
is generic to RME with weak multifractality 77 <gc 1 . 

The level compressibility x m Eq- (H) is obtained by the 
integral of the TLCF 8 ' 10 : 



X=l 



F 2 (e, s) ds. 



(13) 



Using Eq.(0) one can calculate the level compressibility 
X = ^/2 + [l — coth(2/7i)] in the limit of weak multifrac- 
tality 7/ <C 1. Neglecting the exponentially small terms, 
we observe that Eq.(||) (with d = 1) is fulfilled. 

Note that both the linear level number variance £2 (N) 
with x ^ 1 and TLCF of the form Eq.(|J) are not the triv- 
ial consequences of the basis preference. A good counter- 
example is a Gaussian RME with the variance of the 
fluctuating diagonal components different from that of 
the off-diagonal ones. Their ratio /1 = M 2 /X 2 sets the 
new energy scale A > 1 in the problem, such that for 
s 3> A spectral correlations deviate from the Wigner- 
Dyson form. However, the recent analytical results 25,26 
show that this deviation is qualitatively different from 
that described by Eq.(^) for s ^> 1/rj. Albeit the oscil- 
lations in ¥2(5) die out for s 3> A, there is still left a 
constant tail ^(s) = l/7r 2 A 2 that extends up to s = A 2 . 
Therefore the level number variance T,2(N) = (N/nX) 2 
for N < A 2 and T, 2 (N) = N for N > A 2 . 

With increasing 7 and h or decreasing B the fractal 
dimensionality D2 decreases and Eq.(^|) is no longer valid. 
It is reasonable to assume that in the limit h, 7 — > 00 or 
B — > the fractal eigenvector becomes infinitely sparse, 
D p — » 0. For RME-I this is, indeed, the case 14 . In 
this limit Eq.(^) predicts x = 1/2- Let us check this 
prediction using an exact form of TLCF given in Ref. [16]. 

First of all we note that even after unfolding, the TLCF 
R(s, s) in RME-II is not invariant under a shift in s. In 
the limit 7 — > 00 the TLCF has the same form in the 
orthogonal, unitary, and symplectic ensembles 19 ' 22 : 



Y 2 (e, 



-0(1/4 



(14) 



where — 1 < AS < 1 is a deviation of 4e 3> 1 from the 
nearest odd-integer, and 9(x) is a step-function. 

The lack of translational invariance is a peculiarity of 
the particular RME-II. Only the TLCF smoothened by 
averaging over S can be physically meaningful. So we 
arrive at the TLCF of the triangle form: 



Y 2 (s) 



2\s\ 




1, H<l/2, 
otherwise. 



09 = 1,2,4). (15) 



It is remarkable that after the substitution into Eq.(|lJ 
of either Eq.(|l4|) or Eq.(15) we get the predicted value 
X = 1/2. One may thus expect 27 that the triangle form 
of the TLCF is generic for all critical states in the limit 
of an infinitely sparse fractal D2 — > 28 . 
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Let us discuss the applicability of Eqs.(|9]),(|l5|) to the 
description of spectral statistics at the Anderson tran- 
sition. The point is that for RMEs considered here all 
states are critical. For disordered electron systems there 
is a mobility edge E* , and the states are nearly critical 
only in its vicinity where ((E) oc \(E - E*)/E*\- v > L. 
Therefore the RMEs considered here correspond formally 
to ((E) — oo and thus v — oo. Indeed, it has been shown 
8 that for finite v the critical TLCF should have a power- 
law tail R(s) = -c d |s|~( 1+ ~\ where c d ~ \/(ir 2 vd). It 
vanishes in the limit v — > oo in agreement with the expo- 
nential decay of TCLF given by Eq.(||). However, even 
for the realistic case v ~ 1 the power-law tail is small due 
to the additional factor ir 2 d. 

In conclusion, we link together three different random 
matrix ensembles suggested recently. Since in one of 
them the eigenfunction statistics is known to be mul- 
tifractal, we argue that all three RME belong to the 
same universality class with the multifractal eigenfunc- 
tion statistics. By combining known solutions for all 
three ensembles we suggest the form of the two-level cor- 
relation function in the region of weak and extremely 
strong multifractality. 
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